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Deviations from Fermi liquid behavior are well documented in the normal state of the cuprate super- 
conductors, and some of these differences seem to be related to pre-transitional features appearing at 
temperatures above T c . The observation of a pseudogap, e.g. in ARPES experiments, is a familiar 
example of this physics. One potential explanation for this behaviour involves preformed pairs with 
finite lifetimes existing in the normal state above T c . In this way two characteristic temperatures 
can be established. A higher one T* at which pairs begin to form and the actual T c at which a 
phase-coherent superconducting phase is established. 

In order to test these ideas we have investigated the negative U Hubbard model in two dimensions 
in the fully self-consistent ladder approximation at low electron densities. In the non self-consistent 
version of this theory the system always shows an instability towards Bose-condensation of infinite 
lifetime pairs. In contrast to this, pairs obtain a finite lifetime due to pair-pair interaction and the 
sharp two-particle bound state is strongly lifetime broadened when self-consistency is applied. A 
quasi-particle scattering rate which varies linearly with temperature is also found. 
The fully self-consistent calculation we were able to perform using a k-averaged approximation in 
which the self-energy loses its k-dispersion due to a k-average. This approximation is found to 
preserve the essential physics. 
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I. INTRODUCTION: 

In three dimensions (3D) the negative U Hubbard 
model has two well understood borderline cases: One is 
the weak coupling limit U -C t (t is the transfer) . In this 
case the model shows superconductivity of mean-field, 
BCS type with a large coherence length of the cooper 
pairs. The second well understood case is the limit of 
strong coupling and low densities U 3> t, n » 0. In this 
case pairs of electrons which now form composite Bosons 
condense into a two particle bound state and the Fermi 
surface is destroyed. 

In 2D however the situation is different. Here it has 
been shown by Schmitt-Rink et al. ffl that for any cou- 
pling strength the Fermi surface is lost and Bose conden- 
sation (at T=0) takes place. The calculation that was 
used in Q was a T-matrix approximation (taking into 
account ladder diagrams in the particle-particle chan- 
nel) in its non self-consistent, non conserving form which 
is only valid in the zero density limit, n — ► 0. 

The motivation to study a 2D system with low but fi- 
nite densities and intermediate coupling strength (U ~ 
bandwidth W) comes from the high-T c cuprates. In these 
systems the normal state transport is governed by 2D 
Cu02 planes, the quasiparticle density is low but finite 
(one still finds a Fermi surface) and the coherence length 
of a cooper pair is small (£ «3-4 lattice constants) . How- 
ever the negative U Hubbard model in its simplest form 



shows only s-wave pairing. But we believe that full un- 
derstanding of the simple s-wave problem is a necessary 
condition to extend the calculations to more complicated 
models. We are confident that many parts of the physics 
developed here for the s-wave case will survive when more 
complicated d-wave pairing is considered. 

II. CALCULATION 

When trying to expand the known results for the zero 
density limit towards finite densities the main additional 
interaction which has to be taken into account is the 
interaction between pairs. The lowest order interaction 
term would be the exchange of two electrons between two 
pairs. Such interactions are included into the equations 
by extending the non selfconsistent work jO to a fully 
selfconsistent as has been discussed in e.g. 00] • In order 
to perform such a fully selfconsistent calculation we use 
dynamical mean field theory Q which in our particular 
problem becomes not only exact in infinite dimensions 
but also for the limit of large correlations (U > t) which 
is in our particular problem even more important. 

In this case the k independent self-energy (we denote 
the k-average with over-lined quantities) 

= - ^ r(iw m + iuj n )G(iu> m ) (1) 
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results from a K independent vertex function T(iil n ). 
The vertex function now consists of two different parts: 
The one particle continuum and the bound state. In 
the large U limit the K-dispersion of the bound state 
vanishes and therefore the K averaged vertex function 
will already be a good approximation in two dimensions. 
To calculate this vertex function from the susceptibility 
we apply another approximation 



(i - ux(in n )) 
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By doing this the next term which is neglected is of the 
order of the mean square deviation of the susceptibility, 
X 2 — X 2 as is discussed in detail in || . To get the full set 
of equations to solve selfconsistently we further need the 
equation for x 
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and for the one particle Green function: 
G(k,iw„) = (^(k.iwn)- 1 -E(iw„))' 



(3) 
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When solving these equations non selfconsistently we re- 
produce the results of jlj . That means the Fermi surface 
is lost at low temperatures and we get Bose condensation 
into the two-particle bound state. 



III. RESULTS AND DISCUSSION 

When doing a fully selfconsistent calculation the situa- 
tion changes drastically; the infinite lifetime bound state 
gets strongly lifetime broadened and merges with the one 
particle continuum. We further regain a Fermi surface 
H . In Fig.[l] we have plotted the dispersion of the one 
particle density of states A(k,u>) where the self-energy 
was obtained via a fully selfconsistent calculation on a 
quasi-2D system with a constant initial density of states. 
We have chosen several k-points along the (1,1) direction. 
The correlation strongly broadens and renormalizes the 
one particle peak but at the Fermi energy we obtain a 
clear quasiparticle peak. We therefore find no pseudo 
gap from our calculation. 




FIG. 1. For a density of n = 0.3, an attractive correlation of 
U = -8t, and a temperature of ksT = O.lt, the k-dispersion of 
the one particle spectral function is shown. The self-energy 
was obtained by a fully selfconsistent T-matrix calculation 
and the k-points are choosen along the (1,1) direction. Es- 
pecially at the Fermi energy we obtain a quasiparticle peak 
whereas the incoherent broadening comes mainly into play 
away from ki?. The energy units along the axis are chosen in 
units of half the bandwidth [2^1 

From our results it seems to follow that at low but 
finite quasiparticle densities the intuitive picture of pairs 
of electrons which condense as composite Bosons breaks 
down. In the following we discuss some arguments which 
support this result: When thinking in such a Bose picture 
the effective hopping of pairs is given by second order 
perturbation theory. The kinetic energy of such a pair is 
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therefore given by an effective hopping of jj. But only 
in the zero density limit the interaction between such 
Bosons can be considered to be small. At finite densities 
the dominant interaction is given by the fact that a pair 
has a smaller number of virtual hopping processes due to 
the existence of other pairs. Such interaction is therefore 
caused by the Pauli principle, is repulsive can not be 
neglected in comparison with the kinetic energy of the 
pair. Such kind of discussion is well known in nuclear 
physics When one therefore maps the Hamiltonian 
with strongly interacting Fermions onto Bosons one ends 
up, at finite density, with strongly interacting Bosons 
which does not solve the problem. 
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